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ABSTRACT 


Surface combatants face a wide range of threats and perhaps the most destructive 
of these is the underwater explosion generated by a mine or a torpedo. The shock wave 
generated by an underwater explosion can cause severe damage or even a catastrophic 
failure. As the shock wave hits the ship, its energy is transmitted through the structural 
members of the ship. The purpose of this thesis is to examine how this energy is 
transmitted through the watertight bulkhead of a DDG and how the welded stiffeners 
affect the bulkhead’s energy damping properties. To investigate the effects of the 
welding, the bulkhead was modeled both as a finite element model and as a scaled 
physical model. The modes and natural frequencies of the bulkhead were first calculated 
using PATRAN and NASTRAN. Using a one half scale model, the bulkhead was excited 
using random noise over a 250 Hz frequency span while measuring the accelerations at 
60 points along the panel. These measured accelerations were then used to calculate the 
frequency response of the bulkhead and the damping ratios as a function of frequency. 
By plotting the damping ratios versus frequency for each measurement point on the 
bulkhead, there can be a better understanding of how energy waves propagate through a 
welded structure. 
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I. INTRODUCTION 


A. BACKGROUND 

Naval surface combatants face a wide range of threats from the air, the surface 
and below the surface of the ocean. One of the most destructive and perhaps most 
probable threats comes from underwater explosions generated by detonating torpedoes or 
mines. The underwater explosions generated by these weapons create shock waves that 
race out towards a ship at nearly five thousand feet per second. When the shock wave 
reaches the hull of a ship, some of the energy is reflected back into the water but more 
importantly, some of that energy is transmitted to the ship and its structural members. As 
this energy travels through the ship it causes structures to oscillate or vibrate over a range 
frequencies. Vibrations occurring at resonant frequencies can result in large amplitude 
displacements which can cause severe damage or even catastrophic failures which is 
important to understand their effects on a ship’s critical systems and structures. 

Of key concern to the Navy is how a ship and its individual components can 
dissipate the energy imparted by the shock wave of an underwater explosion. This 
dissipation of energy (either over time or distance) is known as damping. The ability of a 
structure to dissipate energy can be affected by many things but primarily by friction at 
structural joints. Construction of Navy ships today relies heavily upon welding, as 
such, understanding the effects welded members have on a structure’s ability to dissipate 
energy is of critical importance. 

In a system such as a steel watertight bulkhead it is clear that the damping is 
related to the material properties but it has also been hypothesized that stress 
concentrations present in the system will have an affect on the overall damping. 
According to Betts, Bishop, and Price, the stress concentrations “account for large local 
increases in the alternating field stress...thereby contributing markedly higher damping 
than the surrounding material that is subject only to the field stress.Additionally they 
hypothesize that fatigue cycling will cause a build up of “exceptionally high local 
damping.In the watertight bulkhead, the stress concentrations are the result of the 
shrinkage that occurs in the area around the welded stiffeners. The shrinkage induces 
tensile yield stress parallel to the weld. The area of these residual stresses has been 
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shown to extend between AVi to 6 plate thicknesses either side of the weld line which can 
account for approximately 15% of the plate material, meaning that the overall damping of 
the panel is increased. 

In a test performed by Betts, Bishop, and Price a mild steel bar, clamped at the 
ends, and held in tension in a Denison machine, was subjected to flexural symmetric 
vibrations (of small amplitude). The damping of the panel was measured by an analysis 
of the decay of free vibration recorded by an accelerometer attached to the center of the 
bar. Plotting the loss coefficient versus mean stress showed that the mean stress appears 
to be a “relevant parameter in deciding the level of damping”^^^ and that “damping is 
significantly increased in the presence of a mean tensile yield stress.” As a result of their 
experiment, the authors believed that “welding effects, together with stress 
concentrations, are among the most important sources of hull damping in deformation 
modes. 

B. OBJECTIVES 

The objectives of this research are to investigate what effect welding has on a 
beam stiffened watertight bulkhead in a DDG. In particular the purposes of this research 
are to determine if the presence of welds in a structure increase the overall damping of 
the structure, does the damping ratio increase at the measurement point as the number of 
welds between the measurement location and the excitation source increases, and to 
compare the Rayleigh damping coefficients from experimental data to those estimated 
from ship shock trial data. 


C. TEST SPECIMEN DESCRIPTION 

This investigation was specifically targeted at the damping characteristics of the 
watertight bulkhead in a Flight lA Arleigh Burke Destroyer. Specimen 1 was a 
rectangular steel panel (102” x 18” x 0.25”) made of A-36 steel, without any welds or 
stiffeners. The unwelded panel was used to gather baseline damping data, how does the 
damping ratio change in relation to distance from the excitation source. Specimen 2 was 
constructed to be a scale model of the watertight bulkhead. Using the structural drawings 
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provided by NAVSEA, a one half scale model (9’ x 4’) made of A-36 steel, was created 
using 9 gauge steel plating, 11 gauge transverse stiffeners, and 12 gauge vertical 
stiffeners. The following tables and figures give the dimensions and characteristics of 
each specimen. 

Specimen 1 - Unwelded Panel 


Table 1.1 Dimensions of Unwelded Panel 


A-36 Steel 

Length 

102 in 

Width 

18 in 

Thickness 

0.25 in 


◄-102” -► 


18 ” 


▼ 


Figure 1.1 Drawing of Unwelded Panel 



Figure 1.2 Picture of Unwelded Panel 
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Specimen 2 - Welded Watertight Bulkhead 


Table 1.2 Dimensions of Watertight Bulkhead 


A-36 Steel 

PLATE 

Length 

9 ft 

Width 

4 ft 

Thickness 

9 Gauge 

IVERTICAL STIFFENERS | 

Length 

47.6875 in 

Width 

2.5 in 

Thickness 

12 Gauge 

IHORIZONTAL STIFFENERS | 

Length 

9 ft 

Width 

2.5 in 

Thickness 

11 Gauge 


9’ 0” 


(” 



_ 13 V4'[ _^ 





Figure 1.3 Drawing o 

f Watertig 

ht Bulkhead 
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II. THEORY 


A. ASSUMPTIONS 

When a system sueh as a ship or a particular component of a ship is excited by an 
external force, it will respond by vibrating and will continue to do so until the force is no 
longer applied. Once the exciting force stops, the system will over time return to its static 
(non-moving) state. The ability of a structure to dissipate the input energy over a period 
of time or across the body of the structure itself is a phenomenon known as damping. 

Total damping in a system can be the combination of many different types of 
damping, including viscous damping, structural damping, and Coulomb damping. 
Although in real structures, many of these forms of damping may be present, a fairly 
accurate mathematical model of the damping can be created assuming only viscous 
damping is present. Assuming viscous damping makes calculating the damping 
properties of the system more convenient equations than if structural or Coulomb 
damping were considered, while at the same time still providing good approximations of 
the true total damping. 

Viscous damping is “the dissipation of energy that occurs when a particle in a 
vibrating system is resisted by a force that has a magnitude proportional to the magnitude 
of the velocity of the particle and direction opposite to the direction of the particle.”^^^ 
Viscous damping can be expressed by the following equation 

F^=cx (2.1) 

c = constant of proportionality 

Another key assumption in creating a mathematical model of the damping that 
can be conveniently calculated is to assume that the system experiences linear 
proportional damping. In other words, the damping of the system is linearly proportional 
to the mass and the stiffness of the system. For a multi-degree-of-freedom system such 
as a ship or watertight bulkhead, these assumptions lead to the following damping 
equation; 

[C] = a[M] + p[K] (2.2) 
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These assumptions form the basis for the model by which the damping values 
across the watertight bulkhead will be calculated. 

B. BASIC DAMPING THEORY 

A vibratory system such as a ship or a watertight bulkhead comprises “means for 
storing potential energy (spring), means for storing kinetic energy (mass or inertia) and 
means by which the energy is gradually lost (damper).The governing equation of 
motion for a damped single degree of freedom (SDOF) system experiencing free 
vibration is: 

mx + cx + kx = F(t) (2.3) 

where 

m = Mass constant 
c = Damping constant 
k = Stiffness constant 


The homogenous solution to the differential equation of motion, Equation (2.1) is 
found by assuming x = e'‘ and F(t) = 0. Substituting into the differential equation yields; 

{^ms^ +cs + k^e'‘ =0 (2.4) 

is satisfied for all values of t when 


2 c k ^ 
s + — s + — = Q 
m m 

Known as the characteristic equation, equation (2.5) has two roots: 


(2.5) 


■^ 1,2 


2m 


V 2m y 


k 

m 


( 2 . 6 ) 


The general solution of the differential equation is thus given by; 

x = AF'‘+BF^‘ (2.7) 

where A and B are constants evaluated from the initial conditions of displacement and 
velocity at time zero; x(0) and x(0). Substituting (2.6) into (2.7) 


X = 


Ae^ ^ + Be ^ ^ 


( 2 . 8 ) 
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The value of the damping term (c/2m)^determines the energy dissipation properties of 
the system. If the damping term is more than (k/m), then the exponents in (2.8) are real 
numbers and no oscillations will occur, a case known as overdamped. If the damping 

term is less than (k/m), the exponent becomes an imaginary number, ±i^jklm-{c ! Imf't , 
and results in oscillations that will diminish exponentially with time, a case known as 
underdamped. For most real structures, the damping ration is almost always less than 10 
percent and is therefore underdamped, thus the mathematical models for this thesis are 
also assumed to the underdamped case. Finally, if the damping term is equal to {k/m) and 
the radical is equal to zero. This case is known as critical damping, Cc 



(2.9) 


As a result any damping can be expressed in terms of the critical damping by a 



C. FREQUENCY RESPONSE FUNCTION 

The purpose of this thesis is to find the damping ratio of the watertight bulkhead 
as a function of time. In order to calculate the damping ratios a modal analysis will be 
performed, but in order to perform the analysis it is first necessary to determine how the 
structure will respond to an input stimulus. For the modal analysis it is more convenient 
to express the equation of motion, equation (2.3) in the frequency rather than time 
domain. The transformation is easily performed using a Fourier transform of the 
equation of motion in the time domain. The resulting equation of motion in the 
frequency domain is given by 

{-mo/+ jco) + k]X{o)) = F{o)) (2.12) 
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or 


where 


and 


X{o)) = H{o})F{o)) 


H{o)) = 


1 

-mo/ + jcco + k 


]_ _ 1 _ 

k \-{a>/a>// +i2(^{a>Icoj 






cOn = Natural Frequency 
= Damping Ratio 
Cc = Critical damping coefficient 


(2.13) 


(2.14) 


(2.15) 


From Equation 2.3 it can be seen that the response of the system { X(o)) }, is 
directly related to the system forcing function {F(o))}, by {H(o))}. If the forcing 
function and the response are known, then { H(o))} known as the Frequency Response 
Function (FRF) can be simply calculated by: 

X(o)) 




F(co) 


(2.16) 


The roots of the characteristic equation, Equation (2.5) are complex roots even 
though the ERE is a function of a real-valued independent variable co. These complex 
roots are also known as the modal frequencies. Rewriting the ERE as a function of the 
modal frequencies gives: 


H(o)) = 


1/m 


(2.17) 


where 

/ =-0),Ci + jO)i^l-Cl 


A partial fraction expansion of the ERE leads to: 

A A* 


H(a>) = 


- + - 


(jco-/) (jo)-/) 


(2.18) 
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The A term and its complex conjugate A* represent the modal constant or mode shapes of 
the structure. 

The above FRF is valid for a system with a single degree of freedom, for more 
complex systems with N degrees of freedom (NDOF) requires the governing equation of 
motion to be described using the mass, stiffness, and damping matrices. 

[M]{v} + [C]{i:} + m{4 = {i^(0} (2.19) 


The individual FRF’s for the NDOF system is developed in the same way as the FRF for 


the SDOF system and the overall FRF is constructed by linear superposition, resulting in: 






+ 




( 2 . 20 ) 


where 

a = Frequency Variable 
p = Measured Degree-of-freedom (Response) 
q = Measured Degree-of-freedom (Input) 
r = Modal Vector Number 
Apqr = Residue 
Xr = System Pole 

n = Number of Modal Frequencies 


D. MODAL ANALYSIS 

Modal analysis is “the process of determining the modal parameters (natural 
frequencies, damping factors, modal vectors, and modal scaling) of a linear, time- 
invariant system.” In order to perform an experimental modal analysis on a structure, 
four basic assumptions are made:^^^ 

1. The structure is assumed to be linear. 

The response of the structure to any combination of forces, simultaneously 
applied, is the sum of the individual responses to each of the forces acting alone. 

2. The structure is time invariant. 

The parameters that are to be determined are constants. 

3. The structure obeys Maxwell’s reciprocity. 
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A force applied at a degree-of-freedom p causes a response at degree-of-freedom 
q that is the same as the response at degree-of-freedom p caused by the same 
force applied at degree-of-freedom q. 

4. The structure is observable. 

The input-output measurements that are made contain enough information to 
generate and adequate behavioral model of the structure. 

The first step in understanding how the modal parameters of a system are 
extracted is to understand the response of a single degree-of-freedom system to a 
disturbance. By calculating the response function for a system, its characteristics or 
modal parameters can be calculated for any given input. 

E. EXTRACTION OF MODAL PARAMETERS 

One method used to extract the modal parameters of the system is the Complex 
Exponential Method (CEM). The CEM requires the use of the time domain version of 
the system response which is found by taking the Inverse Eourier Transform of the ERE 
or the Receptance (displacement/force) «(<») of a viscously damped system. 

n A A* A 

H(m) = aUo) = 2 ,. ( 2 . 21 ) 

The corresponding Impulse Response Function in the time domain is given by: 

2N 

h(t) = Y,A^e^’' ( 2 . 22 ) 

r=l 

Additionally, the Mobility (velocity/force) Y(cl>) of the system can be related to the 
Receptance: 

Y(a>) = ja>a(a>) (2.23) 

The velocity form of the IRE is expressed as: 

2N 

h(t) = Y,AKe"'‘ (2.24) 

r=l 

Erom the velocity form of the IRE, the sampled velocity data set can be expressed as 
follows: 

= h(0),h(At),h(2At),....,h(qAt) (2.25) 

To simplify the following equations, let: 
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( 2 . 26 ) 


. T 7 

e 

Thus the j-th sample data of Equation (2.14) is expressed as: 

2N 2N 

= h(jAt) = Y,AW.j^t) (2.27) 

r=\ r=l 

If extended to the full data set of q sample data, Equation (2.17) becomes; 

A ~ A A A An An 

i\ — VjA A '^AAA ^^2nAnAn 

A = A A + A'A A + • • •+A^ A^ (2-28) 

K = A + V/A A + • • •+ vAAnAn 

Provided that the number of sample points q exceeds 4N, the above equation can 
be used to set up an eigenvalue problem, the solution of which yields the complex natural 
frequencies contained in the parameters Vi, V 2 , etc. 

Multiplying each equation in (2.18) by the coefficient yields the following 
equations 

A) A “ AoA AoA ^ ^ AoAw 

AA = AAA + AAA + • • • + AA^ A^ 

AA = PA" A + A 2 A'A + • • •+ PAInAn (2-29) 

PA = W A + PA2A + • • •+ PAAAn 

Summing the equations in (2.19) results in; 

q IN q 

(2-30) 

(=0 7=1 i =0 

The coefficients Aj taken as the coefficients in the equation; 

Ao + PA + p2y" + Aal^' + • • • + AA" = 0 (2.31) 

With roots Vi, V 2 , Vq. 

By finding the values of the A coefficients, the roots of the system, and hence the 
natural frequencies, can be determined. Recalling that q is the number of data points 
from the IRE and 2N is the number of degrees-of-freedom of the system’s model, by 
setting q=2N; 
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or 

~h\ (2.35) 

L J2Afi2iV ( )2Nx\ 

The unknown coefficients /3 can be found from (2.24) and the roots from equation 
(2.21), thus the natural frequencies can be found using the following relationship; 

y, = (2.36) 

Finally by using equation (2.18) the corresponding modal constants Ai, A 2 ,'", A 2 n 
can be calculated and written as; 



(2.37) 


or 

[V]{A] = {h] (2.38) 

To verify the modal parameters calculated using the above Complex Exponential 
Method, the originally measured time histories are compared to the synthesized or 
regenerated time histories. From equations (2.11), (2.12), and (2.13) 


14 




+ (2.39) 


jCD-X: 


or 


ha,) = t^ 

](0-A^ 




^ , for r > N 


(2.40) 


F. RAYLEIGH DAMPING^^’ 

When analyzing multiple degree-of-freedom systems it is important to consider 
the effect of damping on the complex frequencies and modal vectors. Although one or a 
combination of the aforementioned damping cases may be present in a particular 
structure, it is again more convenient to consider damping that is linearly proportional to 
the system mass and stiffness matrices. This case of damping is known as proportional or 
Rayleigh damping case and has been characterized by the following equation: 

[C] = a[M] + /3[K] (2.41) 

Performing a coordinate transformation to diagonalize the system mass and stiffness 
matrices also diagonalizes the system damping matrix. 

Applying the mass normali z ed modal matrix [^] to equation (2.31) yields; 


Wf [C]W] = = ccl + 

Using equation (2.32) for all 2N modes 

a + = 2(Z)14'i 

a + Pcol = 2 coI (^2 

cc + ~ 



(2.42) 


(2.43) 


(2.44) 


Finally the parameters a and P from equation (2.34) are calculated by matrix pseudo¬ 
inverse; 

{^} = ([«'£.„ ])" [W']L„ (2-«) 
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G. HALF-POWER POINT METHOD^''^ 

Another method for calculating the damping ratio in a system is known as the 
Half-Power Point method. The first step is to record the frequency response function of 
the system in response to an excitation force. For each response peak, the value Q is a 
measure of the sharpness of the peak. When a resonant frequency is reached, cojco^ =1 


and the amplitude peak is given by; 



(2.46) 


The “half-power points” are defined as the points located on either side of resonance 
where v=0.707vmax- Letting v=0.707vmax and squaring the expression for frequency 
response amplitude yields; 


1 


v2^y 


1 


- 

y ^ 


2 

r 


1- 

CO 



2^ 

CO 






v®«y 






Solving for 


f 

' a> ' 


K^nJ 


=(i-2(r*)±2f,/rT 


Assuming □ 1 and neglecting higher order terms; 

= 1±2^ 


f 

' CO ' 




If the roots of equation (2.39) are known as co\ and coi then; 

4^ = 


^2 ^1 ^ 2 




0)2 —co^ 


V y 


(2.47) 


(2.48) 


(2.49) 


(2.50) 


The quantity Q is then defined as; 

=— (2.51) 

2^ 

Thus the damping of the system can be approximated directly from a plot of the 
frequency response system. 
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Figure 2.1 Half-Power Point Method 
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III. EXPERIMENTAL SETUP 


A. EQUIPMENT 

Model 
Model 
Model 
Model 
Model 
Model 
Model 
Model 
Model 

■ Unwelded Steel Panel (102” x 18” x 0.25”) 

■ Vi Scale Model Watertight Bulkhead (9’ x 4’) 

■ National Instruments: Lab VIEW 6.1 Virtual Instruments Software 

■ FORTRAN 

■ Microcal Origin 6.0 

■ MATLAB 

B. EQUIPMENT SETUP 

In order to test the panels in an unrestrained condition, they were suspended from 
a wooden support structure by eight nylon tie-down straps with hooks attached to the 
steel rings welded to the top of the bulkhead while the flat panel was suspended by four 
bungee cords. Once the panel was securely positioned, it was marked with evenly spaced 
numbered points where the accelerometers would be attached. A hole was drilled 
through the panel in order to attach the force transducer and the shaker rod. 

Two methods were used to vibrate the panel. The first method was to excite the 
panel using the PCB modally tuned impulse hammer connected to the HP DSA. The 
second method used the Modal 50A exciter positioned on the lateral excitation stand such 
that when the exciter was connected to the panel the connecting rod was level and 
perpendicular to the panel. The input excitation signal was provided by the HP Digital 
Signal Analyzer connected to the Modal 50A exciter via the MB Dynamics Servo 
Amplifier. The accelerometers were numbered 1 to 5 and placed in order on the panel. 
The accelerometers and the force transducer attached to the shaker rod were then 
connected through the PCB power supply to the National Instruments Virtual Instrument 


PCB Piezotronics Power Supply 
PCB Piezotronics Force Sensor 
PCB Piezotronics Impulse Hammer 
PCB Piezotronics ICP Accelerometers 
MB Dynamics Servo Amplifier 
MB Dynamics Exciter 
Lateral Excitation Stand 
HP Digital Signal Analyzer 
Dell Dimension Desktop Pentium III 


483B07 

208C0I 

086B03 

336C04 

SS250VCF 

Modal 50A 

2050A 

35628 

XPS T800r 
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card installed in the Dell Desktop. The first accelerometer and the force transducer were 
also connected to the two channel HP DSA in order to measure the frequency response 
function of the panel and the signal coherence for a comparison with the results measured 
by the Virtual Instrument. 


Output signal 



Figure 3.1 Equipment Setup for Impulse Hammer Testing 



Accelerometers 


Figure 3.2 Equipment Setup for Exciter Testing 

C. PANEL TESTING 

1. Impulse Hammer Testing 

The initial method of testing was to use the PCB modally tuned impulse hammer 
to excite the panel with an impact, representative of a free vibration system. Connected 
to the HP DSA the hammer provided the input signal while the accelerometers measured 
the system response. Various hammer tips were used to find one that would provide the 
best input signal with the least noise, resulting in a strong coherence. After the input and 
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output signals were amplified the HP DSA analyzed the signals to produce plots of the 
frequency response and coherence. The force-exponential window was used to view the 
plots and 10 stable mean averages were used to generate the frequency response. The 
trigger level for the hammer was set at 0.5 volts. 

2. Exciter Testing 

Instead of using the hammer to simulate a free vibration system, the exciter was 
used to simulate a forced vibration system. The panel as first marked with 60 numbered 
points where the accelerometers would be placed as in Figure 3.3. The exciter is at 
location 10. 



As a first estimate of the damping, the “half-power point” method was used. The 
DSA generated an excitation signal or 1 volt random noise over a frequency range of 0- 
250 Hz. Performing a Fast Fourier Transform over 10 averages, the FRF and the 
coherence were plotted. Where a strong response peak was detected, the magnitude of 
and the frequency at which the peak occurred were recorded. Then using the X and Y 
scroll functions of the DSA a line was created 3 dB below the peak of the response. The 
intersection of this line and the response curve yielded the sideband frequencies. 
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Entering the peak and sideband frequencies into an EXCEL template based on the “half¬ 
power point” method yielded the first estimate of the damping between the stiffeners. 

The second method of measuring the damping of the panel was to measure and 
record the Erequency Response Eunction, coherence, magnitude and phase, and the real 
and imaginary data. The data from the accelerometers and force transducer was recorded 
and processed using the Virtual Instrument 7 Channel DSA. Eor the initial run, the HP 
DSA was setup to generate random noise with a peak value of 1 volt over a frequency 
span of 0-250 Hz. A frequency resolution of 0.05 Hz was desired so the VI was set to 
record 5000 lines of data for each parameter over 10 averages. It was found that a buffer 
size of 20,000 and a scan rate of 100 produced the best coherence when recorded. The 
data was recorded as a text or .TXT file. In addition to generating the input signal, the 
DSA was set to record the magnitude and the ERE as measured by the first accelerometer 
and the force transducer. Using the data gathered from the VI, the .TXT file was then 
input into a EORTRAN 97 program that performed the complex exponential method and 
generated the damping values for a given frequency. The regenerated modal ERE and 
phase were then imported into Microcal Origin 6.0 and plotted with the originally 
measured ERE and phase as recorded by the VI. The extracted damping values and 
corresponding frequencies were then exported into an EXCEL spreadsheet where only 
the frequencies and damping values that corresponded to real modes, identified by 
significant peaks and phase shifts in the regenerated ERE and frequency plots, were 
extracted and used to create the damping versus frequency plot. 

The next step was to calculate the Rayleigh damping coefficients. Using 
MATLAB and inputting the real and imaginary data gathered from the VI the damping 
coefficients a and [5 were calculated. Additionally the program generated a least squares 
curve fit of the damping values. 

When the ERE over the 250 Hz range was examined, it was found that the average 
magnitude of the ERE between 0-20 Hz was much less than the average magnitude above 
20 Hz which resulted in a poorly regenerated ERE plot in the low frequency range. To 
correct this, a separate test was conducted by setting the DSA to generate 1 volt random 
noise over a frequency span of 0-20 Hz. The VI was setup to record 400 lines of data to 
record at 0.05 Hz resolution. The same analysis was performed on the 0-20 Hz data as 
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the 0-250 Hz data. To create a more accurate picture of frequency dependent damping, 
the low frequency and higher frequency damping values were spliced together to cover 
the entire 250 Hz frequency span. For each position, data was recorded and processed for 
both the 0-20 Hz and the 0-250 Hz ranges 

For comparison, the flat, unwelded panel was connected to the exciter and the five 
accelerometers were equally spaced along its length. The same testing parameters as the 
bulkhead were established and the data was recorded in both the 0-20 Hz and 0-250 Hz 
ranges. 


1 2 3 4 5 


Figure 3.4 Accelerometer Positions on Unwelded Panel 

D. FINITE ELEMENT MODELING 

In addition to physical modeling and testing, the watertight bulkhead and the 
unwelded flat panel were also modeled as a finite element model in MSC Nastran/Patran. 
The bulkhead was modeled using 22092 quad 4 shell elements while the flat plat was 
modeled using 29376 quad 4 shell elements. The welds in the bulkhead were not 
modeled; instead the stiffeners and flat panel were merged together by “equivalencing” 
the mesh in Nastran/Patran. The panels were analyzed for their mode shapes at 
frequencies ranging from 0-250 Hz. 
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IV. EXPERIMENTAL RESULTS 


A. IMPACT TESTING 

Although tips of varying hardness were used and many attempts were made to 
produce consistent impacts, the coherence values measured for the impact tests were far 
too low to produce accurate and reliable FRF values. As a result of these inaccurate 
results, impact testing was abandoned early in the experimental process in favor of the 
exciter testing which produced better coherence and could be easily reproduced. 

B. EXCITER TESTING 

I. Half-Power Point Method 

The first estimate of damping values for the panel were found using the “half¬ 
power point” method as previously explained. Using the number 1 accelerometer, the 
force transducer, and the HP DSA to plot the magnitude response of the panel, strong 
resonant peaks were identified and isolated. Once isolated the frequency of the resonant 
peak was found as well as the frequencies of the two sidebands at 3dB below the peak 
magnitude. Using the equation 

2 /„ 

Table (4.1) shows the peak frequencies and the frequencies of the sidebands as 
well as the calculated damping values. An effort was made to use the same resonant peak 
as often as possible, however not every location had the same modal frequency so the 
analysis was performed at the resonant frequencies that had the greatest magnitude. 

Table 4.1a “Half-Power Point” Results (Left Half of Bulkhead) 


f1 = 114.58 1 

1 

f2 = 144.15 
fn = 145.01 

i; = 0.00376 

f1 = 160.25 1 

2 

f2= 160.19 

i; = 0.00045 

fn = 160.33 


f1 = 156.23 1 

3 

f2= 155.6 

i; = 0.00376 I 

fn = 156.78 


f1 = 156.23 1 

Hh 

f2= 155.7 

i; = 0.00312 I 

fn = 156.68 

I 


f1 = 151.75 1 

9 

f2= 151.6 
fn = 151.9 

'C, = 0.00099 

f1 = 130.95 1 

10 

f2= 130.58 

i; = 0.00277 I 

fn = 131.3 

1 Shaker 

f1 = 160.18 1 

11 

f2= 160.09 

i; = 0.00050 

fn = 160.25 

1 

f1 = 169.15 1 

12 

f2= 168.96 

? = 0.00107 

fn = 169.33 



II 

ro 

cn 

17 

f2 = 124.31 

0.00149 

fn = 124.69 


f1 = 135.8 1 

18 

f2 = 135.71 

i;= 0.00050 

fn = 135.85 


CO 

<M 

0 

CD 

II 

M— 

19 

f2 = 160.16 

0.00037 

1 fn = 160.28 



20 


f2= 130.73 i;= 0.00269 
fn= 131.43 


f1 = 131.1 1 

25 

f2= 130.7 
fn = 131.5 

0.00305 

f1 = 170.88 1 

26 

f2= 170.6 

0.00168 

fn = 171.18 


f1 = 130.98 1 

27 

f2= 130.63 

0.00277 

fn = 131.35 


f1 = 169.05 1 

28 

f2= 168.8 

0.00148 

fn = 169.3 
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Table 4.1b “Half-Power Point” Results (Right Half of Bulkhead) 


f1 = 135.55 

33 


f1 = 140.2 

41 


f1 = 140.2 

49 


f1 = 148.93 1 

57 

f2= 135.45 
fn = 135.63 

? = 0.00066 

f2= 139.98 
fn= 140.45 

i;= 0.00168 

f2= 139.9 
fn = 140.49 

? = 0.00209 

f2= 148.76 
fn= 149.08 

?= 0.00105 

f1 = 159.93 

34 

f1 = 151.71 

42 

f1 = 140.3 

50 

f1 = 160.25 1 

58 

f2= 159.83 
fn = 160.01 

? = 0.00059 

f2= 151.54 
fn= 151.86 

£;= 0.00105 

f2= 140.01 
fn = 140.61 

? = 0.00214 

f2= 160.16 
fn= 160.33 

C, = 0.00051 

f1 = 159.85 1 

35 

f1 = 151.58 1 

f2= 151.45 
fn= 151.7 


f1 = 170.83 1 

51 

f1 = 140.25 1 

59 

f2= 159.76 
fn = 159.96 

? = 0.00063 

i; = 0.00082 

f2= 170.55 
fn = 171.08 

0.00154 

f2= 139.98 
fn= 140.53 

?= 0.00196 

f1 = 150.83 1 

36 


f1 = 169.28 1 

44 


f1 = 140 1 

52 


f1 = 140.25 1 

60 

f2= 150.65 

;; = 0.00120 


f2= 169.02 

i;= 0.00157 


f2= 139.7 

i; = 0.00214 


f2= 139.98 

?= 0.00196 

fn = 151.01 



fn= 169.55 



fn = 140.3 



fn= 140.53 



The above table shows the results of the “half-power point” testing. Although the 
damping values vary across the panel, the damping is never greater than 0.4%. 

2. Modal Parameter Extraction 

After completing an initial estimation of the damping using the “half-power 
point” method a more accurate and robust method for calculating the damping was used. 
The National Instruments LabVIEW software was used to perform Fast Fourier 
Transforms on the input data from the force transducer and the output data from the 
accelerometers. The resulting magnitude, phase, real, imaginary, and coherence data was 
exported to FORTRAN via a .TXT file. The FORTRAN code calculated the modal 
parameters for each position using the complex exponential method explained in Chapter 
II. The extracted modal parameters were then imported into Microcal Origin 6.0 where 
the damping values, original and regenerated FRF and phase curves were plotted. 
Comparing the original and regenerated FRF curves, the frequencies and corresponding 
damping values of the most resonant modes (large, narrow banded peaks) were exported 
to an FXCFF spreadsheet. 
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Magnitude 


POSITION 55 
0-20 Hz FRF Curves 



Frequency (Hz) 

*Blue circles indicate true modal responses 

Figure 4.1 Position 55 0-20 Hz FRF Curves 


27 



































POSITION 55 
0-250 Hz FRF Curves 



*Blue circles indicate true modal responses 

Figure 4.2 Position 55 0-250 Hz FRF Curves 

Once imported into the EXCEL spreadsheet the damping values over a range of 
frequencies as well as at a particular frequency over a range of locations could be plotted. 
Einally, the Rayleigh damping coefficients a and (5 and a least squares curve fit of the 
damping at each location was performed using MATLAB and the theory mentioned in 
Chapter 11. 

a. Damping for Horizontal Positions 

The objective of this investigation was to determine the effects welds have 
on structural damping. If welds do result in higher damping, then measuring the damping 
at a point as far from the exciter as possible, in other words the damping should increase 
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as the number of welds between the location of the source and measurement location 
increases. To test this theory, damping measurements were calculated for the locations 
horizontal to the exciter, as shown in the blue box in Figure 4.3. 



Figure 4.4 shows the damping value of each position over the frequency 
range of 0-250 Hz. 


29 






































MEASURED DAMPING FOR HORIZONTAL POSITIONS 
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Figure 4.4 Measured Damping Values for Horizontal Positions 

It is clear that the damping is dependent on the frequency as the damping 
values decrease exponentially as the frequency of the input increases. However it is 
difficult to determine how much the welds contribute to the damping. To determine the 
effects of the welds, the damping values for these positions were calculated at five modal 
frequencies; 16, 59, 76, 181, and 201 Hz. Plotting the damping at each modal frequency 
versus location (inches from the left edge), gave the following results. 
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Table 4.2 


Damping at Modal Frequencies for Horizontal Positions 


POSITION 

DISTANCE 

2 

6.63 

6 

13.25 

10 

19.88 

14 

26.50 

18 

33.13 

22 

39.75 

26 

46.38 

30 

53.00 

34 

59.63 

38 

66.25 

42 

72.88 

46 

79.50 

50 

86.13 

54 

92.75 


DAMPING FOR 

16 Hz 59 Hz 


GIVEN MODAL FREQUENCY 

76 Hz 181Hz 201Hz 


0.00683 


0.00728 

0.00878 

0.00727 

0.00789 

0.00894 

0.00839 

0.00797 

0.00825 

0.00702 

0.00791 

0.00815 


0.00572 

0.00715 

0.00626 

0.00501 

0.00795 


0.00499 

0.00411 

0.00492 


0.00479 

0.00512 

0.00542 

0.00642 


0.00064 

0.00045 

0.00044 

0.00054 

0.00054 

0.00065 

0.00063 

0.00068 

0.00052 

0.00061 

0.00057 

0.00054 

0.00048 


0.00235 

0.00214 

0.00249 

0.00244 


0.00223 

0.00238 


0.00222 

0.00235 

0.00221 


0.00163 

0.00241 


0.00214 

0.00206 

0.00250 

0.00190 

0.00224 

0.00190 

0.00243 
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MODAL DAMPING FOR HORIZONTAL POSITIONS 
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Figure 4.5 Modal Damping for Horizontal Positions 
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From Figure 4.5 it is clear that the damping is not as strongly affected by 
the distance from the source or the number of welds between the measurement and the 
source, but rather the damping is more dependent on the frequency of the excitation. The 
lower frequency excitation (16 Hz and 59 Hz) produced higher damping values with 
more variation than the damping at higher frequencies (76 Hz, 181 Hz, and 201 Hz) 
which had a much lower level of damping with less variation. 

The calculated Rayleigh damping coefficients are listed in Table 4.3 and 
the best fit curve of Rayleigh damping compared to the originally measured damping 
values is shown for each position in the following plots. 

Table 4.3 Rayleigh Damping Coefficients for Horizontal Positions 


Position 

a 

P 

Position 

a 

P 

2 

0.55383 

5.64E-06 

34 

0.76510 

4.90E-06 

6 

0.67804 

4.48E-06 

38 

0.68501 

3.80E-06 

10 

1.13142 

4.60E-06 

42 

0.68214 

3.43E-06 

14 

1.34924 

5.27E-06 

46 

0.46035 

4.26E-06 

18 

0.66848 

5.06E-06 

50 

0.42024 

3.67E-06 

22 

0.76415 

5.04E-06 

54 

0.48189 

4.24E-06 

26 

0.93907 

4.91 E-06 

58 

0.77911 

6.74E-06 

30 

1.56177 

3.42E-06 



The average a value for the horizontal positions was 0.795 with a standard 
deviation of 0.327 and the average P value was 4.63 lE-6 with a standard deviation of 


8.950E-7. 
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Figure 4.6 Best Fit & Rayleigh Damping Curves for Bulkhead Position 2 
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Figure 4.7 Best Fit & Rayleigh Damping Curves for Bulkhead Position 6 
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Figure 4.8 Best Fit & Rayleigh Damping Curves for Bulkhead Position 10 
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Figure 4.9 Best Fit & Rayleigh Damping Curves for Bulkhead Position 14 
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Figure 4.11 Best Fit & Rayleigh Damping Curves for Bulkhead Position 22 
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Figure 4.12 Best Fit & Rayleigh Damping Curves for Bulkhead Position 26 
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Figure 4.13 Best Fit & Rayleigh Damping Curves for Bulkhead Position 30 
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4.14 Best Fit & Rayleigh Damping Curves for Bulkhead Position 34 
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Figure 4.15 Best Fit & Rayleigh Damping Curves for Bulkhead Position 38 
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Figure 4.17 Best Fit & Rayleigh Damping Curves for Bulkhead Position 46 
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4.18 Best Fit & Rayleigh Damping Curves for Bulkhead Position 50 
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Figure 4.19 Best Fit & Rayleigh Damping Curves for Bulkhead Position 54 
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Figure 4.20 Best Fit & Rayleigh Damping Curves for Bulkhead Position 58 


b. Damping for Vertical Positions 

Damping values of the panel were also taken in the vertical direction in 
the same section as the exciter as shown in Figure 4.21 
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Figure 4.22 gives the measured damping values for each position over the 
0-250 Hz range. 
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Figure 4.22 Measured Damping Values for Vertical Positions 

Similar to the trend in damping values for the positions across the panel, 
the vertical damping values in Figure 4.22 show the same exponential decay as the 
frequency increases. Although there are no welds between any of the measurement 
points and the exciter it is still important to check if the location at which the 
measurements were taken has any affect on the damping values. Plotting the damping 
values as a function of location for the modal frequencies (16, 59, 76, 181, and 232 Hz), 
shows that it is more likely that the damping is more likely a function of the frequency 
than the location on the bulkhead. 
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Table 4.4 


Damping at Modal Frequencies for Vertical Positions 


Position 

Distance 

16 Hz 

Damping at Modal Frequency 

59 Hz 76 Hz 181 Hz 

232 Hz 

9 

6 

0.00823 

0.00629 

0.00040 

0.00203 

0.00278 

10 

18 

0.00728 

0.00626 

0.00044 

0.00214 

0.00287 

11 

30 


0.00550 

0.00046 

0.00244 

0.00175 

12 

42 

0.00689 

0.00510 

0.00046 

0.00251 
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Figure 4.23 Modal Damping Values for Vertical Positions 


Figure 4.23 shows for each position the damping value varies from 0.5- 

0.8% for the low frequencies with some variation as the location changes. For the higher 

frequencies the damping value decreases significantly but the damping values for the 

different locations changes fairly little. 

The calculated Rayleigh damping coefficients are listed in Table 4.5 and 

the best fit curve of Rayleigh damping compared to the originally measured damping 

values is shown for each position in the following plots. 
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Table 4.5 


Rayleigh Damping Coefficients for Vertical Positions 


Position 

a 

P 

9 

1.21250 

5.03E-06 

10 

1.13142 

4.60E-06 

11 

0.48038 

5.09E-06 

12 

1.46744 

3.08E-06 


The a coefficients have a mean value of 1.073 with a standard deviation 
of 0.420 while the P coefficients have mean value of 4.450E-6 with a standard deviation 
of 9.398E-7. 






Figure 4.24 Best Fit & Rayleigh Damping Curves for Bulkhead Position 9 
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Figure 4.25 Best Fit & Rayleigh Damping Curves for Bulkhead Position 10 
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Figure 4.26 Best Fit & Rayleigh Damping Curves for Bulkhead Position 11 
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Figure 4.27 Best Fit & Rayleigh Damping Curves for Bulkhead Position 12 
c. Damping for Panel Center 

Damping values were also calculated for a region in the center of the 
panel, as indicated in Figure 4.28. 
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Figure 4.28 Measurement Locations for Damping in Panel Center 
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Figure 4.29 Measured Damping Values for Panel Center 
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Once again the damping is seen to be strongly related to the frequency of 
the excitation. The damping values appear greatest at low frequencies regardless of 
position and appear the least at higher frequencies. Checking the damping values of 
these positions at specific frequencies (16, 25, 76, 82, and 232 Hz) will show the 
relationship between damping, frequency, welds and location. 

Table 4.6 Damping at Modal Frequencies for Panel Center 


Position 

16 Hz 

Damping at Modai Frequencies 

25 Hz 76 Hz 82 Hz 

232 Hz 

22 

0.00789 

0.00251 

0.00065 

0.00098 

0.00243 

23 

0.00555 

0.00241 

0.00060 

0.00149 


26 

0.00894 

0.00185 

0.00063 

0.00055 

0.00248 

27 

0.00703 

0.00209 

0.00064 

0.00173 


30 

0.00839 

0.00247 

0.00068 

0.00124 

0.00265 

31 

0.00783 

0.00205 

0.00050 

0.00106 

0.00254 

34 

0.00797 

0.00188 

0.00052 

0.00047 


35 

0.00865 

0.00077 

0.00055 

0.00079 
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Figure 4.30 Modal Damping Values for Panel Center 
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Figure 4.30 shows again that the damping values for lower frequencies are 
higher and have more variation that the damping at higher frequencies. Figure 4.30 also 
shows that the position vertically or horizontally at which the damping values are 
measured has little impact on the damping. 

The Rayleigh damping coefficients as well as the least fit curves of the 
damping values are shown in the following table and plots. 

Table 4.7 Rayleigh Damping Coefficients for Panel Center 


Position 

a 

P 

22 

0.76415 

5.042E-06 

23 

0.82174 

2.778E-06 

26 

0.93907 

4.913E-06 

27 

0.47627 

3.784E-06 

30 

1.56177 

3.422E-06 

31 

1.31436 

2.902E-06 

34 

0.76510 

4.903E-06 

35 

0.78859 

2.555E-06 


For the center of the panel, the a coefficients have mean value of 0.929 
and a standard deviation of 0.346, the P coefficients have a mean value of 3.787E-6 and a 
standard deviation of 1.038E-6. 
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Figure 4.32 Best Fit & Rayleigh Damping Curves for Bulkhead Position 23 


50 



























Qhwrwa 


p<3atr2fl 




Figure 4.33 Best Fit & Rayleigh Damping Curves for Bulkhead Position 26 
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Figure 4.34 Best Fit & Rayleigh Damping Curves for Bulkhead Position 27 
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Figure 4.35 Best Fit & Rayleigh Damping Curves for Bulkhead Position 30 
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Figure 4.36 Best Fit & Rayleigh Damping Curves for Bulkhead Position 31 
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Figure 4.37 Best Fit & Rayleigh Damping Curves for Bulkhead Position 34 


poatras 



100 1 » 
FREOUEMCV ftifl 


ISO 



Figure 4.38 Best Fit & Rayleigh Damping Curves for Bulkhead Position 35 
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d. Damping for Unwelded Panel 

As a baseline, the damping values in an unstiffened and unwelded panel 
were also calculated in the same manner as the watertight bulkhead. The measurements 
were taken from the positions indicated in Figure 4.12. 


1 2 3 4 5 


Figure 4.39 Measurement Locations for Unwelded Panel 


MEASURED DAMPING FOR UNWELDED PANEL 



Frequency (Hz) 


Figure 4.40 Measured Damping Values for Unwelded Panel 


The damping trend in the unwelded panel follows the same decay trend as 
seen in the watertight bulkhead although the damping values are considerably less than 
those in the watertight bulkhead. It can also be seen that the damping values change very 
little in relation to the position on the panel. 
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Table 4.8 Damping at Modal Frequencies for Unwelded Panel 


Position 

13 Hz 

Damping Vaiues at Modai Frequencies 

27 Hz 66 Hz 93 Hz 178 Hz 196 Hz 

219 Hz 

1 

0.00227 

0.00112 

7.94E-04 

4.14E-04 

0.00105 

0.00128 

5.71 E-04 

2 

0.00228 

0.00119 

8.23E-04 

4.14E-04 

0.00106 

0.0013 

5.61 E-04 

3 

0.00221 

0.00117 

7.90E-04 

4.13E-04 

0.00106 

0.00128 

5.68E-04 

4 

0.00218 

0.00114 

8.08E-04 

4.19E-04 

0.0011 

0.00129 

5.57E-04 

5 

0.00225 

0.00109 

8.13E-04 

3.99E-04 

0.00103 

0.00129 

5.63E-04 


MODAL DAMPING FOR UNWELDED PANEL 


[♦13 Hz m27 Hz Bee Hz A93 Hz 178 Hz >196 Hz ■2191 hz] 



Position 

Figure 4.41 Modal Damping Values for Unwelded Panel 

Just as in the watertight bulkhead, the lower frequencies produced higher 
damping values at each location as well as greater variation in the damping values, but 
still less variation than in the watertight bulkhead. 

The Rayleigh damping coefficients are given in the following table. As 
can be seen the a and P coefficients for the flat panel are considerably less than the 
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damping in the watertight bulkhead which is also evident in the smaller damping values 
shown in Figure 4.40. 

Table 4.9 Rayleigh Damping Coefficients for Unwelded Panel 


Position 

a P 

1 

0.37970 

1.04E-06 

2 

0.36943 

1.07E-06 

3 

0.37198 

1.04E-06 

4 

0.37341 

1.03E-06 

5 

0.36276 

1.02E-06 


For the unwelded panel, the a coefficients have mean value 0.371 of and a 
standard deviation of 0.006. The P coefficients have a mean value of 1.039E-6 and a 
standard deviation of 1.645E-8. The mean values for both the a and P coefficients are 
significantly less than the mean values for any section of the watertight bulkhead which 
suggests that the watertight bulkhead overall has greater damping abilities than the 
unwelded panel. 


f lAT PANEL POSIT 1 




Figure 4.42 Best Fit & Rayleigh Damping Curves for Unwelded Panel 
Position 1 
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Figure 4.43 Best Fit & Rayleigh Damping Curves for Unwelded Panel 
Position 2 
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Figure 4.44 Best Fit & Rayleigh Damping Curves for Unwelded Panel 
Position 3 
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Figure 4.45 Best Fit & Rayleigh Damping Curves for Unwelded Panel 
Position 4 
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Figure 4.46 Best Fit & Rayleigh Damping Curves for Unwelded Panel 
Position 5 
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C. FINITE ELEMENT MODEL 

The finite element model created in MSC Patran/Nastran was used to calculate the 
modal frequencies and mode shapes of both the watertight bulkhead and the unwelded 
panel. 

Table 4.10 Modal Frequencies of Watertight Bulkhead FEM Model 


Mode 

Frequency (Hz) 

Mode 

Frequency (Hz) 

Mode 

Frequency (Hz) 

1 

2.9332 

21 

121.38 

41 

194.26 

2 

20.16 

22 

122.99 

42 

194.51 

3 

32.799 

23 

124.49 

43 

197.85 

4 

42.726 

24 

127.68 

44 

199.6 

5 

5.327 

25 

130.85 

45 

204.86 

6 

61.175 

26 

131.16 

46 

204.88 

7 

65.058 

27 

132.65 

47 

205.07 

8 

68.158 

28 

138.34 

48 

205.43 

9 

76.836 

29 

151.17 

49 

206.15 

10 

87.557 

30 

151.23 

50 

214.44 

11 

87.925 

31 

153.33 

51 

221.85 

12 

88.139 

32 

153.52 

52 

232.51 

13 

89.263 

33 

158.98 

53 

239.11 

14 

97.409 

34 

161.62 

54 

241.16 

15 

100.54 

35 

164.91 

55 

246.39 

16 

108.96 

36 

166.01 

56 

248.51 

17 

113.62 

37 

177.65 

57 

248.57 

18 

114.73 

38 

179.56 

58 

249.3 

19 

115.31 

39 

179.88 



20 

116.91 

40 

190.86 




Table 4.11 Modal Frequencies of Unwelded Panel FEM Model 


Mode 

Frequency (Hz) 

Mode 

Frequency (Hz) 

1 

4.913 

13 

124.250 

2 

13.596 

14 

152.840 

3 

16.821 

15 

158.160 

4 

26.772 

16 

164.400 

5 

34.338 

17 

169.770 

6 

44.436 

18 

178.290 

7 

53.200 

19 

185.940 

8 

66.603 

20 

196.430 

9 

73.988 

21 

201.410 

10 

93.247 

22 

219.410 

11 

97.220 

23 

222.930 

12 

123.380 

24 

243.220 


A comparison of the modal frequencies generated by the FEM model of the 
unwelded panel and the modal frequencies found using the experimental techniques 
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shows that the FEM model frequencies closely match those extracted using the 
experimental techniques, proving the veracity of the FEM model. 

Table 4.12 Comparison of FEM and Experimental Modal Frequencies 


Mode 

FEM Frequency 
(Hz) 

Experimental 
Frequency (Hz) 

1 

4.913 

4.920 

2 

13.596 

13.450 

4 

26.772 

26.550 

6 

44.436 

44.070 

8 

66.603 

66.130 

10 

93.247 

92.710 

15 

158.160 

157.480 

16 

164.400 

164.490 

17 

169.770 

169.670 

18 

178.290 

177.830 

20 

196.430 

195.620 

21 

201.410 

200.550 

22 

219.410 

219.270 

24 

243.220 

242.260 


For the watertight bulkhead many of the highly resonant frequencies that were 
extracted from the experimental techniques were also present in the FEM model; 
however the FEM model also produced many more modes than were seen in the 
experiment which presents the possibility of the FEM model having generated false 
modes. Although the FEM models provide a good picture of the mode shapes and at 
what frequencies they occur, they do not provide information on the damping properties 
of the panels. 
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V. CONCLUSIONS & RECOMMENDATIONS 


A. CONCLUSIONS 

From the experiments performed in the course of this investigation it has become 
evident that the addition of welded stiffeners adds mass and stiffness to an otherwise flat 
panel and increases the Rayleigh damping coefficients, hence the overall damping of the 
structure. It has also been shown that the contribution of these welded stiffeners to 
increases in local damping values is relatively negligible. In other words the presence of 
a single weld or multiple welds between points on a structure did not significantly affect 
the level of damping at any point. 

According to data recovered from the DDG-53 ship shock trials, the Rayleigh 
damping coefficients for ship structural damping were estimated as a=19.2 and P=2.09E- 
6. These experiments however yielded Rayleigh damping coefficients only as high as 
a=0.3797 and P=1.04E-6 for the unwelded panel and as high as a=1.562 and P=6.74E-6 
for the watertight bulkhead. As a result of small increase in Rayleigh damping, 
particularly for the mass (a) coefficient between the unwelded panel and the watertight 
bulkhead and the large difference between the shock trial Rayleigh damping and the 
watertight bulkhead it is unlikely that the welds are a significant source of ship structural 
damping. 

Eurther research into the field of ship structural damping is necessary and 
investigations should be made into the effects of painting, lagging, cable thru-ways, 
hatches, and other components that make up the structure of a ship. 
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